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In efforts to determine phase transitions in the disintegration of highly excited heavy nuclei, a 
popular practice is to parametrise the yields of isotopes as a function of temperature in the form 
Y(z) = z~ T f(z a (T — To)), where Y(z)'s are the measured yields and r, o and To are fitted to the 
yields. Here To would be interpreted as the phase transition temperature. For finite systems such 
as those obtained in nuclear collisions, this parametrisation is only approximate and hence allows 
for extraction of To in more than one way. In this work we look in detail at how values of To 
differ, depending on methods of extraction. It should be mentioned that for finite systems, this 
approximate parametrisation works not only at the critical point, but also for first order phase 
transitions (at least in some models) . Thus the approximate fit is no guarantee that one is seeing a 
critical phenomenon. 

A different but more conventional search for the nuclear phase transition would look for a max- 
imum in the specific heat as a function of temperature T2. In this case T2 is interpreted as the 
phase transition temperature. Ideally To and T2 would coincide. We invesigate this possibility, both 
in theory and from the ISiS data, performing both canonical (T) and microcanonical (e = E* /A) 
calculations. 

Although more than one value of To can be extracted from the approximate parmetrisation, the 
work here points to the best value from among the choices. Several interesting results, seen in 
theoretical calculations, are borne out in experiment. 
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I. INTRODUCTION 



In studies of phase transitions in the disintegration of highly excited heavy nuclei, a popular path for deducing the 
occurrence of a phase transition is to examine the yields of composites. These are readily available from experimental 
data and hence have been the focus of many theoretical studies (jj . The usual practice is || to use a parametrisation 

Y(z) = z- T f(z°(T-T )) (1.1) 

and extract values of r, a and To, which occur in models of critical phenomena The parameters r and a are critical 
exponents and To is the critical temperature. Alternately, in a microcanonical formalism one would write 

Y(z)=z- T f(z°(e-e )) (1.2) 

Here e = E* /A, the excitation energy per nucleon, and eo would be the phase transition energy. Formulae (1.1) and 
(1.2) assume that the thermodynamic limit is reached. In practice, in the nuclear case we have a finite system that 
disintegrates and thus the above parametrisation is only approximate. Hence the values of the parameters can be 
extracted in more than one way and these values may not be the same. We point out that it is largely the Coulomb 
force that leads different methods to give two different values for eo, labeled as e\ and e[. 
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Alternative but perhaps more common tools for studies of phase transition in other fields of physics are measure- 
ments of compressibility, specific heat etc. Experimental data for specific heat were studied in the nuclear case and 
was indeed the cause of great excitement ||. 

We have therefore two distinct ways of trying to deduce a phase transition energy: from the distribution of com- 
posites as the excitation energy is varied (as explained, even here there can be more than one value) or, what may be 
more difficult but achievable, to locate an extremum of the specific heat. We label the excitation energy at which the 
specific heat maximises as ei- It is not obvious that the values of ei,e[ and e2 are close, although from the seminal 
work of Coniglio-Klein jj] on clusterisation this result could be anticipated. 

We have compared both approaches in the nuclear case in the framework of two models. Although the models are 
very different and each has its own strengths and weaknesses, both reveal the following interesting features. If we 
switch off the Coulomb interaction between protons, the deduced phase transition energies, e\,e' x and are close. 
With the inclusion of the Coulomb force, e\ and e\ begin to diverge. For a nucleus of the size of 197 Au, the case 
we study and for which fragmentation data exist, the difference in the e\ and e\ is significant. Further, one of these 
values stays close to the value at which the specific heat maximises and gives a good measure of the phase transition 
energy. 

The two models we use are the Lattice Gas Model (LGM) and a thermodynamic model . The second model 
is close in spirit to the Statistical Multifragmentation Model of Copenhagen ||. We choose to use a microcanonical 
simulation for LGM. So here the primary quantity is the excitation energy e per particle and a temperature can be 
derived afterwards [0]. For the thermodynamic model, we do a canonical calculation so that the temperature is the 
primary parameter and an excitation energy e per particle can be derived afterwards. 

In section 2 we give details of the LGM calculations. Results of the thermodynamic model are presented in section 
3. In section 4 we investigate the ISiS data within those formalisms. We present summary and conclusion in section 



Numerical techniques for microcanonical simulations with LGM have been published |9|. Calculations are done 
for fixed E — Ae where e is the excitation energy per nucleon. This is the primary quantity for simulations. The 
temperature for each simulation can be calculated from T =< 1E^ in > /3. This is discussed in detail in ||. For more 
discussions about the LGM with Coulomb force we refer to |l0|, section II. Bonds due to nuclear forces are taken 
to be -5.33 MeV between unlike particles and between like particles. The LGM has several drawbacks, the most 
noticeable being the lack of quantum effects, which leads to an incorrect caloric curve near T = 0. The LGM has the 
following advantages not shared by several other models. It includes interactions between composites. It incorporates 
the Coulomb interaction in a much more basic fashion (at the nucleonic level) than several other models. This is very 
important for us since this work points to a new effect brought about solely by the Coulomb interaction. Also the 
LGM produces particle-stable composites so that the complicated problem of subsequent particle evaporation 

is circumvented. 

All calculations reported here are for Z = 79 and N = 118 ( 197 Au). At each total energy we compute averages 
after 50,000 simulated events. We use 9 3 (p/po = .27) lattice sites. 

The extraction of parameters from yields (Eq.(1.2)) merits consideration. Discussions of this parametrisation can be 
found in Q where it is used to model a continuous phase transition in an infinite system. As already stated, one does 
not expect the above parametrisation to be exact except in the thermodynamic limit. For very finite systems as is the 
case with disintegrating nuclei formed in very energetic nuclear collisions, the parametrisation is only approximate 
and is by no means a signature of a critical phenomenon but rather that of a phase transition in a finite system, first 
order or otherwise |Q. Theory demands that z of the yields be not small. In the nuclear case it is also not too big 
since the disassembling system itself is very finite. We limit z between 3 and 17, which is similar to most published 
work on the subject. Since the fit is expected to be only approximate, there is more than one prescription for getting 
the "best" parameters. Here we follow the prescription given in Jl3| : 

(1) If Eq. (1.2) were exact, then at e = eo we would have '}2{lnY(z) — InC + rlnz) 2 — 0, as each individual term 
in the sum would be zero. Of course Eq.(1.2) is not exact and thus the sum above will not be zero at any value of e. 
However, this is a valid question to ask. At any given e how well does the distribution fit a power law and what is the 
value of r that gives a best fit to a power law ? At each e we get a "best" r by a least-square fit, i.e., by minimising 
'^(InY(z) — InC + rlnz) 2 with respect to r and C. The "goodness" of fit is given by the smallness of the sum which 
we define to be \ 2 (there are other definitions of x 2 ) 
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II. RESULTS FROM LGM 




(2.1) 
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Here N is the number of terms in the sum. From this step we have a "best" r and a x 2 vs - e. One obvious choice of 
eo is that value of e where \ 2 is minimum (see |jl4j| ). The value of eo deduced using this criterion will be called e! x . 
While this is quite reasonable, it does not use the scaling property z°(e — eo) at this stage at all. That is left upon 
an optimum choice of a afterwards. The more complicated procedure that is followed below is designed to give better 
scaling properties. 

The "best" r vs. e curve will usually have a minimum which we call T m i n - 

(2) Define q = z a (e — eo); f(q) has a maximum for some value of q = q: f max = /(<?)■ For each z the yield Y(z) as 
a function of e has a maximum at some value of e max (z). At this excitation energy Y(z) max = z~ T f max where f ma x 
is a constant independent of z. This allows us to choose values for r and f max using a x 2 test. 

(3) The value of r found above will be higher than T mm . This means if we look for e appropriate for r, two values 
are available from the r vs. e curve (see Fig. 1). The lower value is chosen as the value of eo- The scaling property 
is badly violated by the other choice. The value of eo chosen by this prescription will be labelled by e\. 

(4) Now that we know eo = &\ and e max (z), the excitation at which each z is maximised, we find form a least 
squares fit the value of a from the condition z a (e max — ei) = const, for all z. 

(5) The scaling law can now be tested by plotting Y(z)z T vs. z a (e — ei). Plots for all z should fall on the same 
graph. 
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Fig.l: In this figure we show "best" r and x 2 vs. e (step (1), section II) and also Cy vs. e for the LGM calculation. 
The horizontal line below 2 is a straight line drawn at the value of r from step (2), section II. The locations ei,e2 
and e[ are the values of e for the scaling fit, maximum of specific heat and the minimum of x 2 ■ The left panel is an 
LGM calculation without any Coulomb force, the right panel is with the inclusion of Coulomb. Notice the increase of 
Ae = C2 — ei when Coulomb is included. 



Fig.l depicts graphs (steps (1) to (3)) with and without the inclusion of the Coulomb force. In the graphs we also 
plot the specific heat per nucleon. We refer to the location of the maximum of specific heat as e^. The lessons from 
LGM that we like to emphasize can be learned from Figs. (1) and (2). In Fig. 1 consider first the no-Coulomb case. 
Here e\ is 11.66 MeV, C2 is 12.61 MeV, and e[ (defined by the minimum of x 2 , the prescription of H], see part (1) 
above) is 13.88 MeV. We regard them as close. With Coulomb e\ drops well below e2 and e[; ei and e\ stay close. 

Of the three energies e\, ei and e[, which one marks phase transition better? Without Coulomb, the Y(z) curves 
at the three e-values are quite similar, but with Coulomb, e\ is clearly in the phase co-existence region and is below 
the phase transition energy (please see Fig. 2). Looking at yields at ei and e[, there is no obvious choice in deciding 
which marks the phase transition point better. However, since an extremum in the value of specific heat is a standard 
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signature of phase transition, our preference is with e-i- 




Fig. 2: LGM simulations for Y{z) vs. z at e\,e2 and e[ without Coulomb (top panel) and with Coulomb (bottom 
panel). In the top panel both e\ and e 2 are near the energy, where a maximum in the yield at the high z side has just 
disappeared. Qualitatively, this marks the phase transition point. But in the bottom panel where the Coulomb force is 
included, e\ marks an energy when there is still a large fragment. Thus this is below the phase transition temperature. 
However e 2 still marks the location when the maximum at the high z side has just disappeared. At much larger e 
values (shown arbitrarily at e = 17.6 MeV and e = 14.8 MeV), Y(z) falls much more rapidly with z. 




Fig. 3: For different isotopes z we plot z T Y(z) against z a (e — ei) in the LGM, where the exponent r is the "best" r 
at e\ and extraction of a, e\ is described in the text. By scaling one means that curves for different z 's coalesce into 
one. This is approximately true for scaling around e\ (left panel) but not around e 2 (right panel). Scaling around e[ 
is worse (not shown). 
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In Fig. 3 we show that the scaling law is rather well obeyed around e±. It is very poorly obeyed around e[. An 
interesting plot is a scaling law around e%. This is also shown in Fig. 3. Of course the scaling around e^ is nowhere 
as good as around e\ but it is still better than around e\ (not plotted). 



III. CALCULATIONS WITH A THERMODYNAMIC MODEL 



Details of the thermodynamic model can be found in several places P],pL[7| . The physics assumption is that com- 
posites are formed at an appropriate temperature at a volume larger than normal nuclear volume dictated solely by 
consideration of phase space. Thus the model is close in spirit to the Statistical Multifragmentation Model of Copen- 
hagen H with the simplification that the freeze-out volume is assumed to be independent of the partitions. This allows 
for very quick computation without any Monte-Carlo simulations. The inputs for this calculation are the following. 
Apart from neutrons and protons, experimental binding energies and ground state spins are used for deuteron, triton, 
3 He and 4 He. No excited states are included for these. For mass 5 and higher we use the semi-empirical formula for 
binding energies with volume term, surface tension term, symmetry energy and Coulomb energy. Excited states for 
composites are included in the Fermi-gas approximation. The Coulomb interaction between different composites is 
included in the Wigner-Seitz approximation jq| . 





Excitation energy (MeV/A) 

Fig. 4: Similar to Fig.l except the calculation is with the thermodynamic model (section III), with a freeze- out density 
0.27po- Again note that with the inclusion of the Coulomb force (WC), Ae = ei ~ e\ increases substantially compared 
to no Coulomb (NC) case. 



Since this is a canonical calculation, calculations are done for fixed temperatures. For each temperature, the average 
excitation energy per particle can be calculated. For comparisons with LGM, figures are drawn with energy as the 
abscissa. Table I gives both the temperature and energy for relevant quantities. 

There is no reason to expect results close to the ones calculated using LGM. For example, e\ is lower than e^ in the 
thermodynamic model but higher in the LGM. Nonetheless, Fig. 4 shows that Ae = e2 — e\ increases significantly 
with inclusion of the Coulomb force (from 1.24 MeV to 5.28 MeV). Again e\ does not seem to mark the point of phase 
transition at all (Fig. 5) and e2 is a much better candidate. This is so in spite of the fact scaling is well obeyed with 
respect to ei and not so well with respect to e2 or e[ (Fig. 6). 
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Fig. 5: Similar to Fig. 2 except the calculation is with the thermodynamic model (section III). Again with (WC) 
and without (NC) the Coulomb interaction, e2 continues to be a better mark for the point of phase transition energy. 
Without the Coulomb both e\ and e-i are acceptable. At higher energies (chosen arbitrarily at e—10.8 MeV and 13.^5 
Me V) the drops of yields with z are much faster. 




Fig. 6: Similar to Fig. 3 except that the thermodynamic model is used and scaling is tested around e.\,e.i and e± '. 



IV. 8 GEV/C 7T ON AU DATA 

It is of interest to check if the conclusions reached in the theoretical models are verified in experimental data Jl5],|l6| . 
Several non-trivial issues need to be clarified before this can be attempted. Because the power-law fit is not exact, 
the extraction of r from data or theoretical calculation has some ambiguity. In previous sections we calculated r by 
minimising the quantity defined as x 2 m eq.(2.1) at each e. But one could also minimise 



G 



x 2 



■T\2 



If one is using experimental data, a more standard practice would be to minimise 17 

(Y{z) -Cz- T ) 2 



(4.1) 



(4.2) 



where often the cr's are statistical errors. The difference in the value of r extracted by minimising Eq.(2.1) or Eq.(4.1) 
can be significant or small. In theoretically calculated values of Y(z) the difference is small. But it is not so in the 
experimental values of Y(z). For the experimental data of 8 GeV/c tt~ on gold, the results of using Eqs. (2.1), (4.1) 
and (4.2) are shown in Fig. 7. One can give crude arguments that minimising (2.1) rather than (4.1) means that 
Y(z) , s of higher z's are preferentially fitted. In ifTsf , Eq. (4.2) was chosen. 

In order to compare the experimental results to the theoretical calculations we have taken the experimental yields 
Y(z)'s, ignored all errors and repeated the calculations described in section II. We show in Fig. 7 the results of 
this analysis. Qualitatively, the results are similar to the theoretical results (Figs.l and 4). In those two figures 
the minimum in r occurs at very low excitation energy for the WC case. For the experimental data the minimum, 
if it exists, is also at a low value, below 1.5 MeV. More interesting is the right panel of Fig. 7 where we plot the 
experimental specific heat and find the maximum in specific heat coincides quite well with the minimum of \ 2 (this is 
e[ ) . This agrement is also quite close to the theoretical predictions. The specific heat was extracted by differentiating 
with respect to T, the experimental caloric curve obtained for the same data set by A. Ruangma et al. @. 
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Fig. 7: The left panel shows the values of r extracted from the ISiS data using three different fomulae: (a) by min- 
imising the right hand side of Eq.(2.1) at each e value, (b) by minimising the right hand side of Eq.(4-1) and (c) by 
minimising the right hand side of Eq.(^.2). In the right panel Cy from data and \ 2 (Eq.(2.1)) calculated from data 
are plotted. The minimum of \ 2 and maximum of Cy coincide within experimental uncertainty. 

We also test the scaling around e[ (Fig. 8) although some qualifying comments need to be made about this figure. 
In experiments, the source size as well as the charge of the thermal source, depend upon e. The scaling law, which 
spans e values on either side of e' 1; assumes constant source size as well constant charge. Thus the scaling law can 
not be directly tested without additional corrections renormalising the yields Y(z) to compensate for changes in the 
source size and charge. This was not done here. (However the r values and the values of \ 2 should be insensitive 
to such changes although we do require that for a given e, the source size and charge remain unchanged. This last 
condition is approximately obeyed.) 

We can summarise the results of the comparison with experimental data as follows. In the data the maximum of 
specific heat (e2) and the minimum of x 2 , e[ both are e « AMeV with a value of r about 2.1 and a about 0.53 . This 



7 



is very close to to the results of Elliott et al. |20|] for the same data set using Fisher's droplet model approach and 
Berkenbusch et al. using a percolation model. For percolation, the excitation energy was defined as the critical 
excitation energy in the sense of a second order phase transition. The thermodynamic model, which has no adjustable 
parameters and only a first order phase transition, reproduces trends of the data very well although both e2 and a 
are higher, 6 MeV and 0.96 respectively. 
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z (e-ej) 

Fig. 8: Scaling behavior associated with the ISiS data, around e\. 

Since source sizes change with excitation per particle, we tested the sensitivity of model calculations with regard 
to size using the thermodynamic model. Table II gives the results. Within the range of variation of relevance to the 
experimental data the changes are small, though not negligible. 



V. SUMMARY AND CONCLUSION 



We have discussed three characteristic excitation energies (equivalently temperatures): two of them have their origin 
in Eq.(1.2) (e\= value of eo which gives best scaling and e'^value of eo where deviation from power-law is the least) 
and the third, e2, is obtained from an extremum in specific heat. In an ideal situation, all three would have the same 
value. In theoretical models they are close if the Coulomb force is omitted. The Coulomb force makes a substantial 
splitting between e\ and e'i, indicating the sensitivity of the extracted thermodynamic quantities to this interaction. 
We find that e[ gives a better measure of the phase transition energy and it stays close to ei- In the experimental 
data that we considered e[ and e2 follow this pattern. 

The quantitative results are dependent upon freeze-out densities and the source sizes but not sensitively so. 

Lastly, we have not discussed the order of phase transition but the model calculations in sections 2 and 3 imply a first 
order phase transition. In this interpretation, depending upon the excitation energy, most fragments are emitted while 
inside the co-existence region of the phase diagram (and possibly the spinodal region) and the extracted "critical" 
excitation energy indicates the boiling point. As pointed out earlier, the boiling point is at a similar excitation energy 
as the critical point found by Elliott et al. from an analysis based on the Fisher droplet model and Berkenbusch 
et al. |2|] based on percolation theory. Further discussions of critical phenomen a pfl|,pT| in disassembly of hot nuclei 
as opposed to first order phase transition in the disassembly can be found in [ p"2] , |22| , |23] l 

A first order phase transition is consistent with recent observations by the ISiS collaboration |l6|,^,^5) of a strong 
increase in fragment production probability, a strong decrease in fragment emission time, and the onset of collective 
radial expansion above 4A MeV of excitation energy, which were interpreted as signatures for bulk emission. It is also 
consistent with the flattening of the caloric curve from which the heat capacity was extracted JllJ . 
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TABLE I. The values of the parameters ei,e2 and e' x , corressponding T\,T 2 and T[\ and r's at these ei,e2,ei at ei as 
obtained in LGM and thermodynamic models, for freeze-out density of 0.27po- Values are shown for calculations with (WC) 
and without (NC) Coulomb interactions. 



Parameters 


LGM (NC) 


LGM (WC) 


THDM (NC) 


THDM (WC) 


ei 


11.66 


3.60 


6.77 


0.75 


Ti 


4.46 


1.95 


7.53 


2.94 


r(ci) 


1.964 


1.771 


2.77 


2.36 


e2 


12.61 


8.38 


8.01 


6.03 


T 2 


4.62 


3.38 


7.715 


6.445 


r(e 2 ) 
ei 


1.824 


1.251 


2.44 


1.55 


13.88 


10.5 


7.54 


4.94 


n 


4.88 


3.82 


7.65 


6.05 


r(ci) 


2.544 


1.851 


2.53 


1.42 
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TABLE II. The values of the parameters for two different sizes of the fragmenting source as obtained in the thermodynamic 
model, for a freeze-out density of 0.27po- 



Parameters 


N=118 


N=101 




z= 


Z=DO 


ei 


0.75 


1.01 


Ti 


2.94 


3.46 


r(ei) 


2.36 


2.35 


e2 


6.03 


5.90 


T 2 


6.445 


6.40 


r{e2) 


1.55 


1.49 


ei 


4.94 


4.91 


T{ 


6.05 


6.05 


r(ci) 


1.42 


1.41 
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